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Continuum reduction in large N QCD enables one to extract physical quantities in the 
N — * oo limit of QCD by working in small physics volumes. The computation of chiral 
condensate is an example of such a calculation. 
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1. Introduction 



Oh, QCD with massless quarks does not have any free parameters, 't Hooft 03 pointed 

out that one might consider the size of the gauge group, N, as a parameter. The 
initial hope was to exactly solve QCD in the N — > oo limit and compute physical 
quantities as a power series in 1/N. Although the exact solution of large N QCD is 
' not yet available, significant progress has been in the area of large N phenomenology 

^ . by a careful study of large N counting El. This analysis has extensively shown that 

1/N is a good expansion parameter and that many experimental results can be 
reproduced by just studying one or two leading terms in 1/N. In addition lattice 
studies of pure gauge theories ^ have shown that large N limits of quantities like 
string tension and deconfinement temperature are reached as early as N = 5 and 
N = 6. If one also solves the theory in the N — * oo limit, then one can compute 
physical observables as a power series in 1/N with no other free parameter. 

It seems quite likely that one can solve for the meson spectra in the large N limit 
of QCD using existing numerical techniques and moderate computing power ^ ^^^ 1 
The basic idea behind the numerical solution is the concept of continuum reduction. 
Continuum reduction enables one to work in a finite physical volume and compute 
observables in the infinite volume theory without any finite volume effects. One first 
takes the N — > oo limit at a fixed physical volume and the resulting theory does 
not depend on the physical volume. The numerical study of large N QCD is further 
simplified by the fact that fermions in the fundamental representation are naturally 
quenched. This is the case as long as the theory only has a finite number of fermion 
flavors in the fundamental representation. 
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1.1. Continuum reduction 

Eguchi and Kawai made the observation that the space-time lattice in the limit 
of N — -> oo can be reduced to a single point provided the global U d (l) symmetries 
that multiplies Polyakov loops by U(l) phases are not broken. This symmetry is 
broken for d > 2 and the equivalence of the loop equations used by Eguchi and 
Kawai is ruined But, the arguments of Eguchi and Kawai also hold for an L d 
lattice where L > 1. That is to say, the infinite space-time lattice can be reduced 
to a finite L d lattice in the limit of N — > oo provided the global U d {l) symmetries 
that multiplies Polyakov loops by U(l) phases are not broken on the l d lattice. 
The only parameter in the lattice theory is the 't Hooft gauge coupling b = -h^ ■ 

The location of the transition poin t b c (L) where one of the U d (l) symmetries is 
broken scales properly with L iSLO. One can therefore define a critical size l c in 
the continuum such that the U d (l) symmetries remain unbroken for I > l c . The 
argument of Eguchi and Kawai will hold for all L as long as we keep b < b c (L) and 
there will be no dependence on L. There will be no dependence on the box size I 
as long as I > l c and this is referred to as continuum reduction. Physical results on 
the lattice can be extracted by working on an L d lattice and keeping b just below 
b c (L). Computations should be done on two or three different values of N to study 
the infinite N limit. Once this is done, it is sufficient to work with two or three 
different L values to study the effect of finite lattice spacing and one will be able to 
extract the results for infinite volume QCD in the limit of large N. 

There is one lattice artifact to be taken into account when setting the parameters 
{L,b < b c (L)}. There is a bulk transition on the lattice in the large N limit at 
bs = 0.36 that is associated with the spectrum of the single plaquette (lxl Wilson 
loop). a One has to set b > bs to obtain the proper continuum limit. The scaling of 
the critical coupling b c {L) in perturbation theory up to two loops is given by 
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Lattice study of the phase transtion where one goes from 0c (phase where b > bs 
and all U d {\) are unbroken) to lc (phase where b > bs and one of the U d (l) is 
broken) shows that a value of A = 3.85 ±0.2 in Eqn.Q) fits the data well with b c (L) 
replaced by its tadpole improved value. 



2. Chiral condensate 

QCD with massless quarks has a chiral symmetry that implies {ipif>) = as long as 
one is in a finite physical volume. But the massless limit and the infinite volume do 



a There is no phase transition for finite N but there is a cross-over that gets stronger with increasing 
TV. 
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not commute. In particular, 



lim lim —(ipip(m)} = £ ^ 




(2) 



in QCD at zero temperature indicating that chiral symmetry is spontaneously bro- 
ken. The theory does not depend on the physical volume in the large N limit as long 
as Z > l c . Therefore, the massless limit has to commute with the infinite volume 
limit. How does chiral symmetry break in large N QCD? The answer lies in the 
N — » oo limit. This limit does not commute with the massless limit. In particular, 



and this £ does not depend on I as long as I > l c . 

Overlap fermions ^ enable one to study chiral symemtry breaking away from 
the continuum limit. The computation of the chiral condensate was performed as 
follows. Gauge fields were generated at a fixed L, N and coupling b = . Gauge 
fields were generated in the Q = and Q = 1 topological sectors. Two lowest non- 
zero eigenvalues of the massless hermitian overlap Dirac operator, Ai and A2, were 
computed on each configuration. 13 If chiral symmetry is spontaneously broken as 
N — > 00 at a fixed L and b, then the distributions of z t = Ai£iVL 4 should be given 
by universal functions described by chiral random matrix theory ^5 It was shown 
that the distribution p^zi/zi) was universal as predicted by chiral random matrix 
theory as N gets large at a fixed L and b. Furthermore, we determined £ such that 
both p\{z\) and ^2(^2) followed universal distribution functions. We showed that £ 
did not depend on L after one has taken the large N limit at a fixed b. This confirmed 
that continuum reduction hold for fermionic operators. Furthermore, £(&) obeyed 
the usual scaling law of QCD and was a constant. The detailed analysis 

resulted in the following estimate for the chiral condensate. 



3. 't Hooft model 

It is interesting to consider two dimensional QCD where one has a chiral condensate. 
Continuum reduction is trivial in d = since l c — 0. Therefore, one should be able 
to extract the chiral condensate by working on a l 2 lattice. The only degrees of 
freedom are two SU(N) matrices we will call U and V. The 2N x 2N hermitian 
overlap Dirac matrix, H Q is a dense matrix arising out of the denseness of U and V. 
This model looks very close to a chiral random matrix model with the complication 
arising out of the fact that one has two NxN matrices. We are not aware of an exact 
solution of this model that results in universal distributions of the eigenvalues of H a . 



lim lim _(^(m, N, I)) = S ^ 
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b More eigenvalues can be computed but this increases the computational cost and two are sufficient 
to obtain an estimate of the chiral condensate. 
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Let ±Ai, i = 1, • • • , N be the 27V eigenvalues of H . Then the joint disitribution of 
z.i = x . iN , should obey the universal functions given by chiral RMT since £ = —i — . 

V67TO V67T0 

We performed a numerical simulation of this model for 6 = 1 and several different 
N values. We computed all eigenvalues of H and obtained an estimate of Sj from 
the two lowest eigenvalues by forcing the average of z% to be the same as the one 
dictated by chiral RMT. Fig. ^ shows Ei and £2 as a function of N and one can 
clearly see the approach to the correcy value as N — * 00. Furthermore, the average of 
the ratio, r = j 1 , also approaches the correct value as shown in Fig.^ Fig.[5]shows 
the distribution of pi{z\), P2{z2) and p(r) obtained from the numerical simulation 
and the comparison with the universal distributions given by chiral RMT. The 
distributions approach the expected results as N — > 00. 

4. Conclusions 

Continuum reduction in the large N limit of QCD enables one to compute physical 
observales at infinite volume by performing numerical simulations on finite physical 
volumes. Computation of chiral condensate was used as an example of such a calcu- 
lation. It would be nice to obtain an analytic derivation of the connection between 
chiral RMT and 't Hooft model shown by numerical simulations in section 3. 
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Fig. 1. Extraction of Si and £2 from numerical simulations along with the average value of 
r = i 1 - as a function of N. 



February 1, 2008 4:43 8th workshop on non-perturbative QCD talk 




Fig. 2. Comparison of the distributions of the scaled eigenvalues z\ and 22 with the chiral RMT 
distributions for the same. Also shown is the comparison of p(r) with that of chiral RMT. 



